A Truncated RQ - Iteration forLarge Scale Eigenvalue CalculationsD
نویسندگان
چکیده
We introduce a new Krylov subspace iteration for large scale eigen-value problems that is able to accelerate the convergence through an inexact (iterative) solution to a shift-invert equation. The method also takes full advantage of an exact solution when it is possible to apply a sparse direct method to solve the shift-invert equations. We call this new iteration the Truncated RQ iteration (TRQ). It is based upon a recursion that develops in the leading k columns of the implicitly shifted RQ iteration for dense matrices. Inverse-iteration-like convergence to a partial Schur decomposition occurs in the leading k columns of the updated basis vectors and Hessenberg matrices. The TRQ iteration is competitive with the Rational Krylov Method of Ruhe when the shift-invert equations can be solved directly and with the Jacobi-Davidson Method of Sleijpen and Van der Vorst when these equations are solved inexactly with a preconditioned iterative method. The TRQ iteration is related to both of these but is derived directly from the RQ iteration and thus inherits the convergence properties of that method. Existing RQ deeation strategies may be employed directly in the TRQ iteration.
منابع مشابه
A New Inexact Inverse Subspace Iteration for Generalized Eigenvalue Problems
In this paper, we represent an inexact inverse subspace iteration method for computing a few eigenpairs of the generalized eigenvalue problem Ax = Bx [Q. Ye and P. Zhang, Inexact inverse subspace iteration for generalized eigenvalue problems, Linear Algebra and its Application, 434 (2011) 1697-1715 ]. In particular, the linear convergence property of the inverse subspace iteration is preserved.
متن کاملLarge Scale Empirical Risk Minimization via Truncated Adaptive Newton Method
We consider large scale empirical risk minimization (ERM) problems, where both the problem dimension and variable size is large. In these cases, most second order methods are infeasible due to the high cost in both computing the Hessian over all samples and computing its inverse in high dimensions. In this paper, we propose a novel adaptive sample size second-order method, which reduces the cos...
متن کاملBounds for the Rayleigh Quotient and the Spectrum of Self-Ajoint Operators
The absolute change in the Rayleigh quotient (RQ) is bounded in this paper in terms of the norm of the residual and the change in the vector. If x is an eigenvector of a self-adjoint bounded operator A in a Hilbert space, then the RQ of the vector x, denoted by (x), is an exact eigenvalue of A. In this case, the absolute change of the RQ j(x) (y)j becomes the absolute error for an eigenvalue (x...
متن کاملBounds for the Rayleigh Quotient and the Spectrum of Self-Adjoint Operators
The absolute change in the Rayleigh quotient (RQ) is bounded in this paper in terms of the norm of the residual and the change in the vector. If x is an eigenvector of a self-adjoint bounded operator A in a Hilbert space, then the RQ of the vector x, denoted by (x), is an exact eigenvalue of A. In this case, the absolute change of the RQ j(x) (y)j becomes the absolute error for an eigenvalue (x...
متن کاملFrom O(kN) to O(N): A Fast Complex-Valued Eigenvalue Solver For Large-Scale On-Chip Interconnect Analysis
In general, the optimal computational complexity of Arnoldi iteration is O(kN) for solving a generalized eigenvalue problem, with k being the number of dominant eigenvalues and N the matrix size. In this work, we reduce the computational complexity of the Arnoldi iteration from O(kN) to O(N), thus paving the way for full-wave extraction of very large-scale onchip interconnects, the k of which i...
متن کامل